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Abstract 

We consider the effect of string inhomogeneities on the time dependent background of Brane 
Gas Cosmology. We derive the equations governing the Unear perturbations of the dilaton- 
gravity background in the presence of string matter sources. We focus on long wavelength 
fluctuations and find that there are no instabilities. Thus, the predictions of Brane Gas Cos- 
mology are robust against the introduction of linear perturbations. In particular, we find that 
the stabilization of the extra dimensions (moduli) remains valid in the presence of dilaton and 
string perturbations. 
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I. INTRODUCTION 



Understanding the behavior of strings in a time dependent background has been a 
subject of much interest and has been pursued in a number of differing ways. One scenario, 
known as Brane Gas Cosmology (BGC) is devoted to understanding the effect that string 
and brane gases could have on a dilaton-gravity background in the early Universe ^ O, 
y, y, ^ Q . In P , it was suggested that the energy associated with the winding of strings 
around the compact dimensions would produce a confining potential for the scale factor 
and halt the cosmological expansion^. The analysis of BGC was initially performed at 
the level of the equations for homogeneous and isotropic cosmology. The results were 
recently extended to the case of anisotropic cosmology in j^. There, it was shown that 
string gases can result in three dimensions growing large and isotropic due to string 
annihilation while the other six dimensions remain confined. In jsl it was shown that by 
considering both momentum and winding modes of strings, the six confined dimensions 
can be stabilized at the self-dual radius, where the energy of the string gas is minimal. 
This result demonstrated that, in BGC, the volume moduli of the extra dimensions can 
be stabilized in a natural and intuitive way. 

To date, however, all analyses have been performed at the level of the equations for 
homogeneous cosmology. The string sources are usually represented by a perfect fluid with 
homogeneous energy and pressure densities given by the mass spectrum of the strings (see 
e.g. One may worry that inhomogeneities of string sources (in particular strings 

winding around the confined dimensions) as a function of the unconfined spatial directions 
could lead to serious instabilities which would ruin one of the main successes of BGC, 
namely the prediction that three directions become large, leaving the other six confined 
uniformly as a function of the coordinates of the large spatial sections. 

As a first step towards addressing this concern, we here study the evolution of linear 



^ This was later shown quantitatively in 
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fluctuations about a BGC background. We introduce linear matter and metric fluctua- 
tions, derive the linear perturbation equations, and study their stability on long wave- 
length (wavelengths larger than the Hubble radius) ^ Note that since the dilaton plays a 
vital role in the stabilization process, we must also consider the effect of inhomogeneities 
in the dilaton field. Our main result is that the long wavelength fluctuations do not grow 
in time - there are no instabilities at the level of a linearized analysis. 

In Section 2, we present the equations of BGC and review the background solution 
that leads to 3 + 1 dimensions growing large and 6 dimensions being stabilized at the 
self-dual radius. In Section 3, we introduce the fluctuations and derive the equations 
for linearized perturbations about the time dependent homogeneous background of BGC. 
The long wavelength solutions are given in Section 4, leaving a detailed description for 
the Appendix. We conclude with a discussion of the results in Section 5. 

II. BACKGROUND SOLUTION 

Our starting point is the low energy effective action for the bulk space-time with string 
matter sources j3[, 

S = ^ I d''x^e~'^{R + AiVipr - ^H') + S^, (1) 

where R denotes the Ricci scalar, g is the determinant of the background metric, ip is the 
dilaton fleld, and H is the fleld strength of an antisymmetric tensor fleld. The action of 
the matter sources is denoted by 5*111. For example, with D = 10 this is the low energy 
effective action of type II-A superstring theory. For the purposes of this paper we will 
ignore the effects of branes, since it will be the winding and momentum modes of the 

^ On shorter wavelengths, we expect that the motion of the strings will smear out potential instabilities 
in a way analogous to how the motion of light particles ("free-streaming") leads to a decay of short 
wavelength fluctuations in standard cosmology (see e.g. Q for a review). 
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string that ultimately determine the dimensionality and stability of space-time Here, 
we will ignore the effects of fluxes ^, i.e. we set H = 0. 
This action yields the following equations of motion, 

R + 4V, V> - 4V«</.V> = 0, (2) 

where V is the covariant derivative. 

We will work in the conformal frame with a homogeneous metric of the form 

where (77, x*) are the coordinates of 3 + 1 space-time and x"^ are the coordinates of the 
other six dimensions, both of which can be taken to be isotropic ^J. The scale factors 
a{r]) and 6(77) are given by A = ln(a) and ly = ln{b). 

We consider the effect of the strings on the background through their stress energy 
tensor 

T^^ = diag{p, -pi, -p^), (4) 

where p is the energy density of the strings, pi {i = 1 .. .3) is the pressure in the expanding 
dimensions and pm {iri = 4 ... 9) is the pressure in the small dimensions (because of 
our assuption of isotropy of each subspace, there is only one independent pi and one 
independent Pm)- 

Strings contain winding modes, momentum modes and oscillatory modes. However, 
since the energies of the oscillatory modes are independent of the size of the dimensions, 
and since the winding modes and momentum modes dominate the thermodynamic par- 
tition function at very small and very large radii of the spatial dimensions, we shall here 
neglect the oscillatory modes. In the absence of string interactions, the contributions to 

^ See Ig for inclusion of fluxes in the scenario. 
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the stress tensor coming from the string winding modes and momentum modes (T^ and 
respectively) are separately conserved, 

v%-=o v'^r- = o. (5) 

The conservation equations take the form 

EA:(p-'"^-prn = o, (6) 

i=l 

where the derivatives are with respect to the conformal time 77, and where for the moment 
we consider 9 independent scale factors. 

Expressing Q in terms of the metric Q and the stress tensor (jlj) , we find the following 
system of equations, 

- 3A" - 6u" + 6AV' - 6z/'2 + if" - AV = SyrM-^e^+^A^ 

-A" + 2A'2 + 6AV' + AV = -SnM-^e^+^^pi (8) 
-ly" + 6u'^ + 2AV' + ip'u' = -87rM-2e^+2Ap^ (9) 
-6A" - 12u" - 24AV' - 42i^'2 - 6A'^ - ip''^ + 2ip" + SAV + 12<^V' = 0. (10) 

The explicit forms of the energy density and pressure were given in jsl^, 

p = 3fiN^^'>e-^^~^'' + 3/iM(3)e-^^-6'^ + GfiN^^^e-^^-^" + QpM^^^ e-"^^-^" , (11) 

p. = _^iv(3)e-2A-6. ^ ^M(3)e-^^-6^ (12) 
= _^A^(6)e-3^-5^ + pM^^'^e"^^-^", (13) 

where /i is a constant, iV(3) and M(3) are the numbers of winding and momentum modes 
in the large directions, and A^^^-* and M*^^) in the six small directions. 



^ The equations here are related to Eq. (18) in 5] by the volume factor V — 6-^^+^", e.g. p = tf 
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We are interested in solutions that stabilize the internal dimensions, while allowing 
the three large dimensions to expand. Such solutions were discussed in p, where it was 
shown that the winding and momentum modes of the strings lead naturally to stable 
compactifications of the internal dimensions at the self dual radius, while the other three 
dimensions grow large since the string modes winding them have annihilated, allowing for 
expansion. Thus, we will set A^*-^-* =0. At the self dual radius, the number of winding 
modes is equal to the number of momentum modes (i.e. A^*^^^ = M^^^) and the pressure 
vanishes {prn = 0). 

In Ref. p|, the solutions subject to the above conditions on the winding and momentum 
numbers were found numerically. In this paper, we wish to study the stability of these 
solutions towards linear perturbations in the time interval when the internal dimensions 
have stabilized and the large dimensions give power law expansion. In the following 
section, we will derive the equations for the linear fluctuations. The coefficients in these 
equations depend on the background solution. We will use analytical expressions which 
approximate the numerically obtained solutions of j^. We restrict our initial conditions 
so that the evolution preserves the low energy and small string coupling assumptions 
(gs ~ e^^ < 1). 

A approximation to a typical solution of the equations (fTllTm) is found to be of the form 

A(?7) = A;iln(?7) + Ao or a{r]) = ao'q''^ , 

^{r]) = -k2ln{7]) + ^0, (14) 

where the constants ki, ^2, Aq and ipo depend on the choice of initial conditions. We have 
made use oi u = u' = u" = 0, = 0, = M^^\ = 0. Note that in this limit © 
is trivially satisfied. An example of a solution yielding stabilized dimensions and three 
dimensions growing large corresponds to fci = | and A;2 = |. The numerical solution of 
and the analytical approximation used in this paper are compared in Fig. 1, for the 
above values of the constants ki and ^2. 
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FIG. 1: A comparison between the numerical background solutions obtained in 
line) and the analytical approximation used in this paper (green or dark line). 

III. SCALAR METRIC PERTURBATIONS 

In this section we consider the growth of scalar metric perturbations (see e.g. for a 
comprehensive review of the theory of cosmological perturbations) due to the presence of 



string inhomogeneities. We are interested in the case where the fluctuations depend only 
on the external coordinates and conformal time, not on the coordinates of the internal 
dimensions. For simplicity we work in the generalized longitudinal gauge in which the 
metric perturbations are only in the diagonal metric elements ^. Thus, the metric including 
linear fluctuations is given by 

ds^ = e^^^") ((1 + 2(j))dr]^ - (1 - 2^)5ijdx'dx^) - e^^^''^ (1 - 205m„(ix"^rfx". (15) 

The dilaton (/? also fluctuates about its background value lpq. The dilaton fluctuation x 
is determined by 

if = ipo + 5ip X = Sv- (16) 

In the above, the fluctuating fields X) 0) "0 and ^ are functions of the external coordinates 
and time, i.e. 

X = x{v,x'), = 0(r/,x*), ip = ip{T],x'), ^ = ^{r],x'). (17) 

The perturbations of the matter energy momentum tensor result from over-densities 
and under-densities in the number of strings. From ()ll|) - (fT^ and noting that we are 
interested in the case when A^(3) = and M^^^ = N'-^^ we find, 

5p = 5p^ + 6p^, (18) 

Sp^ = 30pN^e-^^ + ISpN^lje-^^ + e^A^^^^e'^^, (19) 

5p^ = 42pN^e-^^ + ISpM^e-^^ + ISpNipe^'^^ + UfiMi/je-'^^ + QpSM^^^e"^^ + (20) 

+3/i5Me-^\ (21) 

5px = QfiM^e-^^ + ApMije-^^ + pSMe-^\ (22) 

Sp^ = 2pN^e-^^ - p5N^^'>e-'^^ + ^i5M^^h-^^, (23) 

^ As discussed e.g. in jlo| . for scalar perturbations depending on all spatial coordinates it would be 
inconsistent to choose the perturbed metric completely diagonal, and one would have to add a metric 
coefficient to the dtdx"^ terms, where x™ are the coordinates of the internal dimensions. However, as 
discussed in Q, if the fluctuations are independent of the coordinates as in our case, the coefficient 
can be chosen to vanish, and thus the perturbed metric is completely diagonal. 
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where we define^ N = A^^^) = M^^) and M = M^^\ The fluctuations 5N^^\ 5M^^\ and 
5M are taken as functions of both conformal time and the external space, e.g. 5N^^'^ = 
5N{ri,x'). 

It follows from (jH)) that the perturbed sources obey modified conservation equations 
for both the winding and momentum modes, 

5p'^^^ + y: \'^{5p-'^ - 5pjn + y: sKip-^-^ - pjn = o, (24) 
1=1 1=1 

where 6X = a^^6a = —ij) and du = b^^6b = — ^ are spatial variations. 

We rewrite (0) to take the more familiar form of the Einstein and dilaton equations, 
namely 

g^'d.^d,^ - \g^'V,V.ip = \e'^Tli, (25) 

where we invoke Planckian units (i.e. SnMp"^ = 1). Plugging the perturbed metric (fT3j) 
and dilaton into these equations, making use of the background equations of motion, and 
linearizing the equations about the background (i.e. keeping only terms linear in the 
fluctuations) yields the following set of equations: 

- - V^x + 37^x' + 20(/^'' - 2xV', (26) 

a,a,(0-V'-6e-2x) =0 1^3, (28) 
^ Notice that we must be careful to distinguish between the perturbed quantities SN^^^^ and 6M^^\ 
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[df - v^) (0 - - 6^) - 2ij" - 6^" - AHij' - en^' - 27^20 - m'(f) - 2(p'n 

+2x" - 2V\ + 27ix! + 40<^'' - 4xV', (29) 

- + 5V^e - + 2v2v^ - SV'" - 107^^' - 30'7^ - 97^v^' - m^4> - m'cj) 

= e2^+2A(^2xr™ + 57-) - 40y." - 20V' " - H'^' - m'^' + 2x" 

-2V2x + 4Hx' + 40<^'2 - 4x'0', (30) 

+ \^\ - W = \e^^+'\2xT + 5T) , (31) 
where T = T^ is the trace. The modified conservation equations take the form 

|^(5iV(6)) =7iVe', (32) 

-12/z(5MA'e-^ + 3/i^ (5M) e"^ = 0. (33) 

These equations give us the evolution of the metric perturbations 0, ■0, and ^ in terms of 
the matter perturbations x, 6p, and At first glance, it may appear that the above 
system is over-determined since we have eight equations for seven unknowns. However, as 
is the case in standard cosmology, the conservation equations are not independent of the 
Einstein equations. Thus, we can choose to keep only one of the modified conservation 
equations and our system will be consistent. 

IV. EVOLUTION OF FLUCTUATIONS 

In this section we focus on the evolution of long wavelength perturbations, since insta- 
bilities on these scales would be the most dangerous for the success of BGC at providing 
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a mechanism which allows exactly three spatial dimensions to become large and stabilizes 
the radius of the other dimensions at a microscopic value. We focus on the time interval 
rj ^ 1 when the hierarchy in scales between the large and the small spatial dimensions has 
already developed, and where we can use the analytical approximations to the background 
dynamics given in Section II. 

Since we are focusing on long wavelength fluctuations, we can neglect all terms with 
spatial gradients in Equations Making use of (PHj) we can eliminate one of 

the scalar metric perturbations, 0, from the rest of the Equations ()2fi|) - (j!?T|) . We take 
Equation as a constraint on the initial data, noticing that it only contains first 
derivatives of the perturbation variables. By using the background solution ()14|1 we find 
that the remaining equations represent a coupled system of harmonic oscillators with time 
dependent coefficients. With these approximations, the equations become 

-2x" - 2r - ee - f - g^' - Igx' -i§^X-^i^-^C = 0, (34) 
-Sr - - 2x" - W,^' - - igx' -^^-m^i-m?X = ^, (35) 

-k" - igx' + l,^' -i^X-£^i^-M^^ = 0. (36) 

Notice that the string source terms, 6N and 6M, have disappeared from the equations. 
This is due to the fact that, for long wavelength fluctuations, the terms due to purely 
gravitational dynamics (e.g. the damping terms in the evolution of the radius of the 
internal dimensions as a consequence of the expansion of the large dimensions) are more 
important than the matter sources. This is not too surprising based on the results of 
the theory of linear cosmological perturbations which show that on super-Hubble scale 
matter is pulled along by gravity, but that it is the purely gravitational dynamics which 
determines the growth rate of the fluctuations. In analogy, we find here that the terms 
representing the string matter sources are sub-leading in the equations (see Appendix). 
This realization is the key physical reason which leads to our ultimate conclusion that 



''' The full equations can be found in the Appendix. 
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also in BGC there are no instabilities of long wavelength fluctuations. 

The above set of coupled equations can be solved exactly, and we find the leading terms 
to be 

^-^3 + ^, (38) 
X ~ C5 + (39) 

This is our main result. First, we have found that the perturbation of the internal ge- 
ometry ^ is constant to leading order. Thus, the stabilization mechanism resulting from 
considering winding and momentum modes remains viable in the presence of string in- 
homogeneities and dilaton fluctuations. Secondly, the geometry of the 3 + 1 dimensions 
growing large does not suffer from exponential instabilities. In fact, we have power-law 
decay for and the dilaton perturbation x decays to a constant. 

The crucial point that led to this behavior was the fact that the string source terms were 
sub-leading in the equations above. Physically this means that the number of winding and 
momentum modes moving in or out of a region is negligible compared to the expansion 
rate and the evolution of the perturbations. Therefore, there is an averaging of the long 
wavelength perturbations that leads to a smoothing effect and the string inhomogeneities 
are rendered harmless. 

V. CONCLUSIONS 

We have considered the evolution of linearized string inhomogeneities and dilaton fluc- 
tuations in the background of Brane Gas Cosmology (BGC). We considered fluctuations 
which are independent of the spatial coordinates of the internal dimensions. We have 
derived the perturbation equations for BGC in the longitudinal gauge and solved them in 
the long wavelength approximation. We have found that for long wavelengths, the effects 
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of the string sources is subleading compared to the purely gravitational and dilaton terms 
in the fluctuation equations. This generalizes the result of the conventional theory of 
cosmological perturbations in Einstein gravity in four space-time dimensions which states 
that on scales larger than the Hubble radius, self-gravitational effects dominate the dy- 
namics, and that the matter sources are simply dragged along by the metric. In particular, 
there are no instabihties of the background solution of BGC towards such long wavelength 
fluctuations. Thus, we flnd that the predictions of BGC are robust towards the effects of 
matter and dilaton fluctuations. In particular, the stabilization mechanism for the extra 
dimensions (volume modulus) remains operative in the presence of these inhomogeneities, 
while the 3 + 1 dimensional space-time continues to grow without instabilities. 

Although this is a promising result for BGC, there are still many questions to be 
answered. Our analysis has been limited to the linear regime, while non-linear effects 
may be crucial. Our linearized analysis corresponds to taking the localized string sources 
and smearing them out over a length scale larger than the Hubble radius to obtain a 
string gas description. It would be interesting to see if these results hold when considering 
individual strings and their interactions and in what limits the string gas approximation 
holds. Also, BGC offers a dynamical way to stabilize the volume modulus, but we have 
not addressed the issue of stabilizing the dilaton and we have also ignored the role of 
fluxes. We leave these questions for further investigations. 
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APPENDIX A: PERTURBATION EQUATIONS 

In this Appendix we present a detailed derivation of the solutions presented in section 
4. In the long wavelength limit we can neglect gradient terms and the equations pH|l - (j!?T|) 
reduce to 

-mi' - sHtp' - sn^ = |e2^+2^(2xro + 6t^) - mcj)^' - sv^v' - eev + 

+3W + 20¥^'2 - 2xV, (Al) 
a,a,(0-^-6e-2x) =0 1^3, (A2) 

-2^Ij" - Qi" - mij' - QHC - 2n^4) - m.'4> - 20'7^ = e^^+'^'p{2xTi + 5Tl) + 
-40(^" - 2(/)V' - ^'H4>ip' - AtP'^' - 12^V' + 2x" + 27^x' + ^^'^ - 4xV', (A3) 

-5f' - 3< - lOHi' - - QHij' - en^cl) - QH'(t) = e^^+^^ (2xT^ + 6T^) , 
-4(j)ip" - 20V' - 8ncl)ip' - 6^V' - lo^V' + 2x" + mx' + - 4xV', (A4) 

-20^'2 + + W - G^W + - 27^0¥^' - I^V - seV - \x" - nx' 

^ ig2A+2^(^2x/ir; + 5r;). (A5) 

Making use of the background solution (fT^ and ()A2|) to eliminate the scalar metric 
perturbation (/), we find the following set of equations describing the evolution; 

61 ^/ 61 f / 88 „// I ( 11114 3000 12000 _L I 36000 9000 11114 \ t j_ 

-2T^X -7^4 -21^^ + y-T323n^-;^-^jX+y-^-^-m^jt + 

+ (-^ - W - ™ V - -^^M('^ - -^SN('^ - -^SM(') = 0, (A6) 

\ 77"BT id^.ir? J 2»7"BT 

9 v" - 97/)" fi;^" 230 f / 176 // _ 230 / I ( _ 2000 _ 12868 ) ^ i ( _ 4000 _ 6434 ] / , 
-^X - 0^ - - iMiW 21^X + I ^ 3969;^ / V / 

+ - iwV - ^^^^'^ = 0, (AT) 
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5344 



V 

5344 



jX + -^SN^^^ - -^6M('^ = 0, (A8) 

V 



1 ^ 11 250 , / , 9 J J , 170 6000 , 9000 43 K/, i 510 15000 1 c i 

~ 63^^ + 7^^^ + I ^49^^ ~ 1^ + V ^ " 1^ + I ~"4V ^ + 



3 -5Ar(6) _ 0. (A9) 



The first equation, containing only first derivatives, is taken as a constraint on the 
initial data. We are interested in the late time behavior, i.e. 77 ^ 1. Keeping only the 
leading order terms the remaining equations can be approximated as, 

-2x" - 2r - er - - w,^' - w,x' -i^x- - = 0, (mo) 

-3r - Sr - 2x" - ig^' - Ife' - if X' - lig.^ - - Tig.X = 0, (All) 

iv" - 250^/ , A,// - J^v 43_ ; _510_t _ n ( AT9'\ 

This set of coupled equations can be solved exactly and we find the leading terms to 

be 

^-^3 + ;^, (AM) 
X-c, + ^^. (A15) 

As a consistency check, we can plug this result back into the original equations ()A6j) - 
flA9|l . In addition we must consider the evolution of ^A^'-^-', SM^^'^ and 5M given by (jH^ 
and ()33|). We find that it was consistent to neglect the time evolution of the string sources 
(i.e. dN and 6M) compared with the expansion of the background and the evolution of 
the perturbations. In addition, we find that at late times this holds as an exact solution of 
the perturbation equations. In this way, we find a self consistent solution for the behavior 
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of the long wavelength linear perturbations. 
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